Medium induced jet absorption in relativistic heavy ion collisions 



Axel Drees, 1 Haidong Feng, 1 and Jiangyong Jia 1 ' 2 

department of Physics and Astronomy, State University of New York, Stony Brook, NY 1 1794-3800 
Current Address: Columbia University, New York, NY 10027 and Nevis Laboratories, Irvington, NY 10533, USA 

(Dated: February 8, 2008) 

The dense medium created in Au + Au collisions at the Relativistic Heavy-Ion Collider (RHIC) 
significantly suppresses particle production from hard scattering processes and their characteristic 
back-to-back angular correlation. We present a simple model of jet absorption in dense matter which 
incorporates a realistic nuclear geometry. Our calculations are performed at the jet level and assume 
independent jet fragmentation in the vacuum. This model describes quantitatively the centrality 
dependence of the observed suppression of the high pr hadron yield and of the back-to-back angular 
correlations. The azimuthal anisotropy of high pr particle production can not be accounted for 
using a realistic nuclear geometry. 

PACS numbers: 27.75.-q 



I. INTRODUCTION 

In relativistic heavy ion collisions at 62.4 < ^/snn < 
200 GeV in the Relativistic Heavy-Ion Collider (RHIC), 
partons scattered with large Q 2 become a valuable tool 
to study nuclear matter under extreme conditions. Due 
to the large Q 2 , the hard scattering occurs early in the 
collisions, thus the scattered partons may directly probe 
the subsequently produced hot, dense and strongly inter- 
acting medium. 

In collisions of elementary particles, i.e. in the absence 
of a dense medium, the hard scattered partons typically 
fragment into two back-to-back jets of hadrons. These 
hadrons have large transverse momentum (j>t) an d pro- 
nounced azimuthal angular correlations. The production 
cross sections can be calculated with perturbative QCD 
(pQCD) based on the fragmentation theorem and extrap- 
olated to heavy ion collisions. Theoretical studies predict 
that the hard scattered partons suffer a large energy loss 
by multiple scattering and induced gluon radiation as 
they propagate through dense matter [l| . Unlike energy 
loss in QED, the radiative gluon energy loss per unit 
length dE/dx not only depends on the color charge den- 
sity and the momentum distribution of the partons in 
the medium, but also linearly depends on the thickness 
of the medium, due to the non-Abelian nature of gluon 
radiation in QCD HQ. 

Data from Au + Au collisions at ^/snn = 130 and 200 
GeV have revealed rich information on high px phenom- 
ena. One of the most interesting results is the apparent 
"jet quenching" in central collisions, observed as suppres- 
sion of the hadron yield by a factor of 4-5 [J 0. @ and 
hadron back-to-back correlation strength by a factor of 
5-10 0, compared to expectations based on the under- 
lying nuclcon- nucleoli collisions. The absence of these 
phenomena in d + Au collisions suggests that the ob- 
served suppression in central Au + Au collisions is in- 
deed an effect of the dense medium created during the 
collisions,consistent with parton energy loss in the dense 
medium &B1 EES- 

Both the suppression of the high px hadron yield and 



the back-to-back angular correlations show a characteris- 
tic centrality dependence, which seems to be independent 
of pr for pr > 4.5 GeV/c Interestingly, in this 

Pt range, particle production seems consistent with jet 
fragmentation, despite the suppression. Specifically, ex- 
periments have observed: 

1 . An identical spectral shape compared to p + p col- 
lisions within systematic errors |^, ||| . 

2. A similar suppression for charged hadrons and 
7r 0, s [f| and for charged hadron, A, and [T^ . 

3. A h/ir° ratio consistent with values observed in p 
+ p collisions, indicating a similar particle compo- 
sition in p + p and Au + Au at high px 0, El • 

4. Strength, width, and charge composition of near 
angle correlation consistent with jet fragmenta- 
tion 0. 

5. xt scaling of pion production cross section in Au + 
Au between ^Jsnn. = 130 and 200 GeV is similar 
to p + p collisions [j| . 

The data suggest that high px particles come dominantly 
from jet fragmentation in vacuum. 

In relativistic heavy-ion collisions, energy is deposited 
in the overlap region between two colliding nuclei. The 
size, shape and energy density of the medium formed in 
this region strongly depends on the impact parameter of 
the collision. The amount of medium a hard scattered 
parton traverses, and subsequently its energy loss, varies 
with the centrality of the collision and also the azimuthal 
angle with respect to the reaction plane. If the parton 
energy loss is large, the surviving partons will be emitted 
dominantly near the outer layer of the overlap region |14| . 
The partons moving towards the surface (near side) tra- 
verse on average less material than those going in oppo- 
site direction (away side). Thus partons scattered to the 
near side are likely to escape with little energy loss, while 
the away side partons are likely to lose significant energy 
and thus are suppressed more strongly. 
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It was shown before [TEL Il6l Il7l ITsf that for a static 
non-Abelian partonic medium, the mean total energy loss 
depends quadratically on the medium size. In addition, 
the 'dynamical scaling law' allows to write the parton en- 
ergy loss for an expanding medium in terms of an equiva- 
lent static medium with a linear dependence on the path 
length 17]. Thus to the first order, one can decouple 
details of the energy loss formulation, which is compli- 
cated and model dependent, from the simple geometry 
of the medium. Recently, several authors have modelled 
the energy loss of hard scattered partons in Au + Au col- 
lisions InA 

in mm mm. 

Different theoretical ap- 
proaches for the energy loss have been used and compared 
to experimental data like the nuclear modification fac- 
tor Raa, the back-to-back jet correlation strength, and 
the azimuthal anisotropy of particle production V2- In 
general, reasonable agreement between the model calcu- 
lations and experimental data can be achieved. 

In this paper, we present a simple model of jet ab- 
sorption in an extremely opaque medium, following the 
approach in Ref. (2^ that was used to set a limit on the 
expected azimuthal anisotropy due to jet quenching. The 
model is equivalent to an extreme jet energy loss scenario 
in which most of the energy is lost in a single scattering. 
We use this simplified model to focus on the geometric 
aspect of the jet propagation in dense medium. After in- 
corporating a realistic nuclear geometry, we demonstrate 
that the centrality dependence of the high pt suppression 
of the yield and of the back-to-back correlation strength 
can be described quantitatively. We also discuss the sen- 
sitivity of our results on the jet absorption patterns and 
different collision geometry assumptions, and make pre- 
dictions for the system size dependence. 



II. COLLISION GEOMETRY AND JET 
ABSORPTION: THE MODEL 

Our discussion is limited to centrality dependence of 
the suppression pattern over the px range of 4.5 to 10 
GeV/c, where the hadron suppression is approximately 
constant Throughout the discussion, the average sup- 
pression values for each centrality bin from the data are 
used to compare with our calculation. Furthermore, since 
hadron production in this pt region seems to be con- 
sistent with vacuum jet fragmentation, we neglect the 
medium modification of jet fragmentation and assume 
that the suppression for partons is identical to the one 
observed for hadrons 1 . 



Due to the steeply falling power-law spectra of scattered partons, 
the hadron spectra is dominated by fragments biased towards 
large z. As was shown in Ref. 1261 . the shape of the hadron spec- 
tra becomes nearly identical to the shape of the parton spectra, 
independent of the fragmentation process. Since we model jet 
absorption at the parton level, the suppression of hadrons thus 
should be approximately the same as the suppression of partons. 



The collision geometry is modelled by a Monte Carlo 
simulation of Au + Au collisions based on the Glauber 
approach [23]. For the Au nucleus a Woods-Saxon den- 
sity distribution with radius R = 6.38 fm and dif- 
fusivity a = 0.53 fm is used [2^. We calculate for 
each simulated Au + Au collision the underlying num- 
ber of nucleon-nucleon collisions (A^qJj) and the number 
of participating nucleons (-^Vpart)' assuming the inelastic 

nucleon-nucleon cross section is crjy^ = 42 mb. 

Centrality classes are defined from the fractional cross 
section by cutting on the impact parameter of the colli- 
sions. However, the centrality determination is insensi- 
tive to the specific cuts and consequently methods em- 
ployed elsewhere [2^, H(],E1 give similar results. The av- 
erage number of participants, nucleon-nucleon collisions 
and impact parameter for different centrality classes are 
listed in Table. HJ For each centrality class, we also calcu- 
late the participant density profile, Ppart^' v)> an d t ne 
nucleon-nucleon collision density profile, p co ii (x, y), in 
the plane transverse to the beam direction. The density 
profiles for a central event class and a peripheral class 
are shown in Fig. ^ the peak values of the collision and 
participant densities are also listed in Table [I] for all cen- 
trality classes. 

In the following we assume that the energy density is 
proportional to the participant density. This is moti- 
vated by the recognition that the bulk particle produc- 
tion scales approximately with the number of partici- 
pants [3(], [32, l33l l34| . In a later section we study the 
sensitivity of our results on this assumption. To give the 
participant density a physical scale we relate it to the 
energy density using the Bjorken estimate |35l ]. 

E ^part m 

where to = 0.2 - 1 fm/c is the formation time, ro = 1.2 fm 
is the effective nucleon radius, and Ap ar ^ is the number of 
participating nucleons. For central collisions the exper- 
imentally determined value of e^j is 5 GeV/fm 3 j3|| 
for To = 1 fm/c. With approximately 350 participants 
the scale factor to convert participant density to ei-,; is 2 
GeV/fm for r = 1 fm/c or 10 GeV/fm for t = 0.2 fm/c. 

Binary scaling of hard scattering assumes that the in- 
coming parton distribution in Au + Au collisions is a 
superposition of the individual nucleon parton distribu- 
tion functions. According to the factorization theorem 
the probability for the hard scattering processes in Au + 
Au is then proportional to p co \\- Therefore we generate 
back-to-back parton pairs in the transverse plane with 
a distribution of p co y[(x,y) and isotropically in azimuth. 
These partons are then propagated through the nuclear 
medium with density proportional to Ppart^: y)- The 
survival probability of a parton produced at (x, y) along 
direction (n x ,n y ) is calculated as |25| 

/ = exp(-K/) , (2) 
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FIG. 1: (Color online) Collision density p co rj(a;,y) (a) and participant density PpartC^J/) O 3 ) m transverse plane for central 
(0-5%) and peripheral (75-80%) collisions. The bottom panels are the average density along the impact parameter axis, which 
is indicated by the dashed line in the top panels. 



TABLE I: . Centrality dependence of parameters calculated with a Monte Carlo based Glauber approach. For each centrality 
class, we list the average number of participants, the average number of nucleon-nucleon collisions, the average impact parameter, 
the average participant density in the overlap region, the maximum participant density, the maximum nucleon-nucleon collision 
density, the maximum energy density calculated using the Bjorken estimate at to — 1 fm/c(see text), and the transverse area 
with a density larger than 1 GeV/fm 3 . 
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350 
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3.8 
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67 
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35 
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1.09 


1.5 


2.7 


3.1 
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75-80% 


11 


10 


13 


0.79 


0.76 


0.9 


1.5 




6.9 


90-95% 


4.1 


2.8 


14.5 


0.56 


0.3 


0.23 


0.6 








where k is the absorption strength, which is the only free 
parameter in the model. / is the matter integral along 
the path of the parton, which is calculated as, 

I= dl l-m-p( x + ( l + l o) n x,y+(l + lo)n y ) (3) 

JO i + t0 

This parameterization corresponds to a quadratic depen- 
dence of the absorption (oc Idl) in a longitudinally ex- 
panding medium (ri^ - ) 2 • Here we assume that partons 
travel with speed of the light and that they sense the 
dense matter after a formation time of 0.2 fm/c or a dis- 
tance of l = 0.2 fm. However, the results presented in 
this paper do not depend strongly on the choice of Iq. 

The absorption strength n is then fixed in central Au 
+ Au collisions to reproduce the observed hadron sup- 



2 The effec t due to transverse expansion was shown to be 
small |2(3 . 



pression. In the following, we use k = 0.7 which gives 
a suppression factor of 4.35 as measured for 0-5% most 
central collisions by PHENIX Q. This corresponds to 
an absorption length of A ~ 3.4 fm for a parton travers- 
ing an expanding medium with a constant participant 
density of 2.48 fm -2 . 

Fig. Elshows the origination position (x, y) of the hard- 
scattering partons that escape the medium in 0-5% most 
central collisions. The depletion at the center of the over- 
lap region is evident. The origins of the surviving jets 
are almost uniformly distributed over the collision re- 
gion, which biases jet emission towards the surface of the 
collision region. Partons emitted from the center of the 
collision region typically traverse matter with density of 
more than one participant per fm 2 for a distance of > 5 
fm and thus arc frequently absorbed, while partons gener- 
ated near the surface can easily escape if emitted towards 
the surface. Therefore, the short absorption length (less 
than half of the radius of the nuclear overlap) naturally 
leads to emission of jets and thus of high pt hadrons from 
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FIG. 2: a) Origination point distribution in the transverse 
plane for partons escaping the dense medium created in cen- 
tral collisions with (6) ~ 2.2 fm. b) Generated (solid line) 
and survived (dotted line) parton initial x position for the 
cut \y\ < 1 fm. 



outer layer of the medium. 



III. CENTRALITY DEPENDENCE OF JET 
ABSORPTION 



Once k has been determined, the survival rate of jets, 
the probability to find back-to-back jets, and also the az- 
imuthal anisotropy of emitted jets relative to the reaction 
plane is fixed for any centrality class. In Fig. [21 we com- 
pare our model calculation to the centrality dependence 
of three compilations of experimental data: 

1. Charged hadron and vr° data from PHENIX 0,0 
and charged hadron data from STAR Q are plotted 
on the top panels of Fig. [3 Shown are the ratios of 
the observed high px hadron yields relative to the 
expected yield from the underlying nucleon-nucleon 
collisions, normalized to -W^ji, Raa in panel a), or 

ATp ar ^/2, Rj^ V ^ in panel b), of the specific cen- 
trality class. The PHENIX data have a lower px 
cut of 4.5 GeV/c, while the STAR data are given 
for px > 6 GeV/c. 

2. The bottom left panel (c) gives the back-to-back 
jet correlation strength measured by STAR 0. 
The away side correlation strength can be defined 
byH 



D A A(PT ig ) 



(4) 



trig 



dp T 



/'o 



da h A \/cPp^ ig 



for an associated hadron /12 with px in backward az- 

tritj 

imuthal direction of a trigger hadron hi with p T . 
The STAR data are for charged trigger particles 

with 4 < p^ lg < 6 GeV/c and associated charged 
hadrons with px > po — 2 GeV/c detected within 
|0i — 02 1 > 0o = 2.24. The data are normalized to 
the expectation from p + p collisions, and corrected 
for combinatorial random background and the az- 
imuthal anisotropy of bulk particle production. 

3. The azimuthal anisotropy of high px charged parti- 
cles, quantified by V2, the second Fourier coefficient 
of the dN/d(j) distribution, is shown in the bottom 
right panel (d). The PHENIX data are measured 
for 4 < px < 5 GeV/c with respect to the re- 
action plane |37j . The preliminary V2 values from 
STAR experiment were determined using reaction 
plane method around 6 GeV/c and 4 particle 
cumulants for 5 < pr < 7 GeV/c [39| . 

The results of our model calculation are shown as thick 
solid lines on all panels. Once the absorption strength is 
adjusted to the observed yield from the most central col- 
lisions, the centrality dependence of the normalized yield 

^part 

is well reproduced. In Fig. yp, the calculated R AA 
increases for peripheral collision, which is expected if the 
yield scales with the A^jj (thin solid line) . Only a small 
fraction of the partons is absorbed since the matter den- 
sity and volume are small. As the centrality increases 
both matter density and volume increase; for collisions 
with more than 100 participants absorption overwhelms 

the increase due to point like scaling and the R AA 
decreases with centrality. 

The jet absorption model also reproduces the magni- 
tude and centrality dependence of the back-to-back corre- 
lations without further adjustments. The calculated sup- 
pression is negligible for peripheral collisions, increases 
continuously and reaches a suppression factor of 7 for the 
most central bin, consistent with the data within errors 3 . 
The model suggests that the suppression of back-to-back 
correlations is almost a factor of 1.6 stronger than that 
for the single inclusive hadron yield. In the jet absorption 
model, most of the surviving partons come from near the 
surface of the overlap region. Thus their partner partons 
have to traverse on average a longer distance through the 



3 We notice, however, due to the large errors of the STAR data, 
the suppress factor for central collisions is not well constrained. 
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medium. This naturally leads to a stronger suppression 
of the back-to-back correlation. It should be noted that 
in our model, we assume that the jet fragments outside 
of the medium, in other words, the fragmentation is iden- 
tical to p + p by construction. Therefore, the near angle 
jet correlation strength can not deviate from unity in our 
model. 

Due to the asymmetry of the overlap region of the two 
nuclei, the average amount of matter traversed by a par- 
ton depends on its azimuthal direction with respect to the 
reaction plane, which leads to an azimuthal anisotropy of 
the emitted jets This anisotropy reaches its maxi- 

mum for collisions with an impact parameter of about 9 
fm corresponding to approximately 100 participants. It 
is small for peripheral and central collisions. Our cal- 
culation reproduces the measured trend of the centrality 
dependence of v-i, but the magnitude is below the mea- 
sured value 4 . With a matter density profile deduced 
from a Woods-Saxon distribution, a large fraction of the 
surviving jets is emitted from the low density region at 
large radii (see Fig. Eh)- Therefore the anisotropy is 
diluted. Larger values of V2 are obtained for larger re, 
but these reproduce neither the suppression of the yield 
nor the back-to-back correlation strength. For a Woods- 
Saxon nuclear profile, v-i reaches a maximum at a certain 
re (for b = 9 fm, i; max = 10%), but then decrease to 
as re — > oo because all surviving jets come from the dif- 
fuse outer layer of the overlap region, which is essentially 
isotropic [2J| . The calculated v% values are very sensitive 
to the actual nuclear profile used, we shall come back to 
this in Section II V Bl 



h = I dl lp(x + ln x ,y + ln y ) (5) 
Jo 

I2 assumes absorption proportional to the path length I 
in a longitudinally expanding medium as one may expect 
for a hadronic energy loss |4(| scenario; I3 assumes a 
static source and absorption oc I 2 . Again the absorption 
strength is adjusted to give a factor of 4.35 suppression 
of the yield for central collisions. We find re = 0.84 fm 
(A = 1.9 fm) and re = 0.06 (A = 3.6 fm) for I 2 and I 3 , 
respectively. 

The results are compared with each other in Fig. ^ 
The centrality dependence of the normalized yield is 
similar for the three different absorption patterns and 
thus gives little discrimination power. The correlation 
strength D aa is more sensitive, but the largest sensi- 
tivity is found for the anisotropy parameter V2- The I2 
scenario tends to localize the absorption in the region 
within the distance of ~ Iq from the jet creation point. 
In this case, the suppression is dominated by the initial 
matter profile and is insensitive to the later evolution of 
the system. This naturally leads to a similar value be- 
tween Raa and Daa- The surviving jets are also emitted 
more isotropically resulting in a smaller v 2 5 . 

On the other hand, the absorption in a static medium 
with quadratic path length dependence (I3) has a strong 
dependence on the jet path. Partons are absorbed along 
their full trajectory in the medium. Thus both the sup- 
pression of the Daa and the value of i>2 depend more 
strongly on the centrality of the collision. In particular, 
one observes a stronger centrality dependence of Daa 
than of the Raa and a larger V2 ■ Although the magni- 
tude of V2 is still below the experimental value. 



IV. DISCUSSIONS 

A. Dependence on Absorption Pattern 

Parton energy loss through gluon bremsstrahlung is 
thought to be proportional to the path length squared 
in a static medium 0, ■ This motivated our ansatz in 
Eq. [3J We note that although / calculated from Eq. [3| 
was interpreted as absorption oc Z in a longitudinally 
expanding medium it can also be treated as absorption 
oclina static medium when Iq is small. In the following, 
we denote it as I\ . In order to test the sensitivity of the 
data to discriminate between different types of energy 
loss we repeat the calculation for two additional types of 
absorption, 

f°° '0 

h = I dl p(x + (I + l )n x ,y + (I + l )n y ) 

Jo I + <o 



B. Dependence on Density Profile 

Motivated by the approximate iVp ar ^ scaling of bulk 
particle production, we have assumed that the initial en- 
ergy density created is proportional to the participant 
density (e oc Pp^t)- A closer look at the data reveals 
that both the transverse energy Et and particle multi- 
plicity increase faster than linear with -A^p^ |29l l33l I ll| . 
In a two component model 01 1 the additional increase 
is attributed to (mini)jet production, which should scale 
with -A'cqIP To test the sensitivity of our model to the 
initial energy density profile, we repeated the calculation 
assuming that the energy density is proportional to the 
collision density p co \\{x, y). 

The results of the calculation are presented in Fig. 
Generally the agreements between the data and the 
model calculation are equally good for both matter den- 
sity profiles. In detail the suppression of the inclusive 



4 We note that the 6% V2 from our calculations corresponds to a 
factor of 1.3 larger jet survival probability in reaction plane than 
out of plane. 



J In this analysis, we have ignored the transverse expansion, 
which can further reduce V2, but only slightly change Raa and 
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FIG. 3: (Color online) Centrality dependence of normalized yield of high pt hadrons ( a and b ), back-to-back correlation 
strength (c), and V2 (d) at high pr- In all four panels, the thick solid line indicates the result of our calculations based on 
Eas. 12131 Further details are discussed in the text. 



yield (Fig. [SJi and b) shows a slightly stronger central- 
ity dependence. This is easily explained: the absorption 
strength k is fixed to reproduce the suppression of the 
yield in the most central collisions; (p co \\{x, y)) varies 

by factor of 4 faster than (pparti^' 2/)^ fr° m peripheral 
to central collisions. As a result, the absorption is smaller 
in peripheral collisions. Similarly Daa is less suppressed 
in peripheral collisions, but decrease faster for interme- 
diate centralities. On the other hand, the centrality de- 
pendence and magnitude of v-i are similar to Fig. in- 
dicating that the two density profiles produce an almost 
identical anisotropy for all three absorption patterns. 

In some recent model calculations nuclear distributions 
different from a Woods-Saxon profile have been used, typ- 
ically either a hard-sphere [23, [25| or cylindrical nuclear 
distributions 0, 0, 0|. The results are significantly 
different for these unrealistic approximations. In par- 

,. , . . i , cylindrical sphere 

ticular, V2 is increased, where v 2 " > v 2 > 

. .Woods-Saxon 

u 2 

Fig.Elshows the calculations for a hard-sphere nuclear 
distribution in the top panels and for a cylindrical nuclear 
distribution in the bottom panels. In both cases, all three 
absorption scenarios miss the centrality dependence of 



the normalized inclusive yield. In contrast to the data, 
the suppression sets in at rather peripheral collisions and 
remain approximately constant with centrality. 

On the other hand, cylindrical and spherical profiles 
apparently result in much better agreement with vi data. 
The reason is that both cylindrical and hard-sphere col- 
lision profiles have sharp surfaces with a large eccentric- 
ity, leading to a larger u 2 0. ^ n contrast, a Woods- 
Saxon collision profile has a diffuse surface. While it 
may be conceivable that the density profile of the matter 
produced in the collision deviates from the convolution 
of two Woods-Saxon nuclear distributions, it is hard to 
imagine that the probability for hard scattering deviates 
from a Woods-Saxon density distribution significantly. 

C. Energy and System Size Dependence 

If jet suppression is really dominated by final state en- 
ergy loss, it should decrease if the energy density or the 
volume of the medium is reduced. In addition to vary- 
ing the impact parameter, both volume and density can 
be reduced by varying the mass number of the colliding 
nuclei. Alternatively the density may also be changed by 
varying the beam energy ^/snn. 
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FIG. 4: (Color online) Calculation of Raa (a) R AA ^ Xt (b), 
Daa (c), and V2 (d) using a Woods-Saxon nuclear profile for 
Ii (solid line), I2 (dotted line), and 73 (dashed line) types of 
jet absorption. 
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Np ar i dependence of the Raa and Daa is remarkable. 
This implies that to the first order the volume is propor- 
tional to iVp ar t and that the participant densities profiles 
are very similar. This similarity is illustrated in Fig. [S] 
by comparing the integral distribution 7i experienced by 
the generated partons in central collisions for various col- 
lision systems with that for centrality selected Au + Au 
collisions with similar iVp ar t. Since the jet correlation 
strength also depends on the volume and average den- 
sity, the calculated Daa is also similar for light systems 
and corresponding Au + Au collisions. In contrast, the 
anisotropy calculated for central collisions of smaller nu- 
clei is consistent with zero as expected. 

To guide our estimate of the beam energy dependence 
we assume that the matter density scales like the Bjorken 
energy density. For an average mix of quark and gluon 
jets which we do not vary with y/s, we fix the absorption 
coefficient k to be the same as in ^/snn — 200 GeV. In 
this estimation the energy density increases by a factor of 
2 from SPS (^/sWn = 17 GeV) to RHIC (^/sWn = 200 
GeV) energies. Interpolating between these values and 
scaling down the density in our calculation accordingly, 
we find that the high px hadron suppression is still a 
factor ~3 for lower energy RHIC run at y / s/vjv = 62.4 
GeV. 

Our estimation would also predict a factor of 2 sup- 
pression at SPS energies, which is not observed experi- 
mentally. However, data at the SPS ^(| are limited topx 
below 4 GeV/c, a region where other mechanism like pt 
broadening and hydrodynamic flow complicate the inter- 
pretation of particle production in heavy ion collisions. 
We noted that the predictive power of our approach is 
limited, because we did not take into account the ^/snn 
dependence of the quark (gluon)jet fraction, the hard- 
scattering cross section and the system lifetime, which 
could be important. 



CONCLUSION 



FIG. 5: (Color online) Calculation of R AA (a), R AA (b), 
jet correlation (c), and V2 (d) assuming e ~ P co n- The three 
curves are for 7i (solid line), I2 (dotted line), and ^3 (dashed 
line) types of jet absorption, respectively. 



The jet absorption model can provide a simple base- 
line prediction for both the beam energy and system size 
dependence of the suppression, if we assume that the ab- 
sorption strength k is the same as in y / Sjv"/v=200 GeV 
Au + Au collisions. We limit the discussion to the cal- 
culation presented in section HTl 

The system size dependence can be evaluated without 
further assumptions. In Fig. the nuclear modification 
factor Raa for central collisions of smaller systems (see 
Table. |Hj at v /sjv r /v=200 GeV is compared to the calcu- 
lated centrality dependence for Au + Au collisions (see 
Fig. I3J . The agreement of the absolute values and the 



We have shown that the experimentally observed cen- 
trality dependence of the suppression of the hadron yield 
Raa and the suppression of back-to-back correlation 
Daa can be quantitatively described by jet absorption 
in an opaque medium when a realistic nuclear geometry 
is included. In our model, the centrality dependence of 
Raa and of the Daa are rather insensitive to the details 
of the medium density profile. Raa is also insensitive to 
the absorption patterns. Both Raa and Daa do not dis- 
tinguish between central collisions of small systems and 
centrality selected Au + Au collisions with similar Np ar ^. 
Our interpretation is that Raa and Daa depend mostly 
on volume and average density of the opaque medium. 
In the limit of very opaque matter, Raa probably will 
not reveal further insight into the details of the absorp- 
tion pattern or alternatively of the energy loss and the 
density profile, unless data at much higher px become 
available. 
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cylindrical nuclear profile (bottom panels) using the 7i (solid line), I2 (dotted line), and I3 (dashed line) form of jet absorption. 



TABLE II: List of light collisions systems and corresponding Ap ar t for 0-5% central collisions. 



Species 


O + O 


Si + Si 


Fe + Fc 


Cu + Cu 


Zr + Zr 


I + I 


La + La 


^part 


24 


45 


95 


109 


159 


224 


246 




part part 



FIG. 7: (Color online) The calculated system size dependence 
(open circles) of the Raa (left) and Daa (right) compared 
with their centrality dependence in Au + Au collisions (thick 
solid line). The open circles are calculated for 0-5% most 
central collisions for each collision system. The thick solid 
line represent the calculation for 7i type of absorption (same 
calculation shown in Fig. and c.) 



We find that observables like Daa and V2 are more 
sensitive to the absorption patterns. Our model does not 
describe V2 quantitatively unless an unrealistic nuclear 
density profile is used. This might indicate that the ac- 
tual diffuseness of the opaque medium is smaller than 
expected from the density profile obtained by convolut- 
ing two Woods-Saxon nuclear distributions [l!j. It is also 
conceivable that the real suppression is larger, and that 
the observed suppression is reduced by "soft" particles 
from dynamic mechanisms different from jet fragmenta- 
tion, such as hydrodynami cs j48l ] plus viscosity correc- 
tion J43, quark coalescence [47|. and quark/diquark pick 
up |24|. In this case, both the "soft" particles and a 
stronger suppression would lead to a larger V2- 
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